Finding is as easy as detecting for quantum walks 



Hari Krovi 1 , Frederic Magniez 2 , Maris Ozols 3,4 , and Jeremie Roland 3 

1 University of Connecticut 
2 LRI, Univ Paris-Sud, CNRS 
3 NEC Laboratories America, Inc. 
4 University of Waterloo and Institute for Quantum Computing 



Abstract. We solve an open problem by constructing quantum walks 
that not only detect but also find marked vertices in a graph. The number 
of steps of the quantum walk is quadratically smaller than the classical 
hitting time of any reversible random walk P on the graph. 
Our approach is new, simpler and more general than previous ones. We 
introduce a notion of interpolation between the walk P and the absorbing 
walk P', whose marked states are absorbing. Then our quantum walk is 
simply the quantum analogue of the interpolation. Contrary to previous 
approaches, our results remain valid when the random walk P is not 
state-transitive, and in the presence of multiple marked vertices. 
As a consequence we make a progress on an open problem related to the 
spatial search on the 2D-grid. 



1 Introduction 

Many classical randomized algorithms rely heavily on random walks or Markov 
chains. The notion of hitting time is intimately related to the problem of spatial 
search, where displacement constraints are modeled by an undirected graph G. 
The set of desired vertices, or marked vertices, is denoted by M . Classically, a 
simple algorithm to find a marked vertex is to repeatedly apply some random 
walk P on G until one of the marked vertices is reached. The hitting time of P, 
HT(P, M), is precisely the expected number of repetitions, or steps, necessary 
to reach a marked vertex, starting from the stationary distribution of P. 

Quantum walks are natural generalizations of classical random walks. Am- 
bainis [1] was the first to solve a natural problem - the "element distinctness 
problem" - using a quantum walk. Following this, many quantum walk algo- 
rithms were discovered [2-4]. Quantum walk algorithms for the spatial search 
problem [5] were studied for the hypercube [6] and the grid [7,8]. 

The notion of hitting time has been carried over to the quantum case in [8- 
13]. Usually, the quantum hitting time has a quadratic improvement over the 
classical one. However, until the present paper, several serious restrictions were 
imposed: a quantum algorithm could only detect the presence of marked ver- 
tices [10], but for being able to find them the Markov chain has to be reversible, 
state-transitive, and with a unique marked vertex [14,13]. The detection al- 
gorithm is quite intuitive and well understood, whereas the finding algorithm 



requires an elaborate proof whose intuition is not clear. This is due in part to a 
modification of the quantum walk, so that the resulting walk is not a quantum 
analogue of a Markov chain anymore. 

Whether this quadratic speed-up for finding a marked vertex also holds for 
any reversible Markov chain and for multiple marked vertices was an open ques- 
tion. In this paper, we answer this question in the positive. Here we choose 
another approach by modifying directly P, and by considering the quantum 
analogue of the modified random walk. Doing that we keep some intuition and 
get simpler proofs while obtaining more general results. The new walk is sim- 
ply an interpolation between the walk P and the absorbing walk P', whose 
marked states are absorbing. The interpolation coefficient can be used to tune 
the overlap between the stationary superposition of the quantum walk and its 
projection onto marked vertices. For a suitable value, depending on the reparti- 
tion of the marked vertices, the overlap with marked and unmarked vertices is 
balanced, leading to a quantum walk algorithm that finds a marked vertex within 
y/HT(P, M) steps (Theorem 5). The balancing can also be achieved when lim- 
ited or even no information is available on the repartition of the marked vertices 
(Theorems 6, 7 and 8). As a consequence, we make a progress on an open 
problem from [5, 14] related to the spatial search on the 2D-grid (Corollary 2). 

2 Preliminaries 

2.1 Spatial search on graphs 

We fix throughout the paper an undirected graph G = (X, E) . Let M C X be a 
set of marked vertices and let m = \M\. Vertices are encoded in a distinguished 
vertex register. Our goal is to find any of the marked vertices in M using only 
evolutions that preserve the locality of G on the vertex register, i.e., to perform 
a spatial search on G [5]. Here we define an even more restricted notion of 
locality than the ones in [5], but it is more intuitive and sufficiently powerful 
for our purpose. We allow two types of operations on the vertex register: static 
transformations (that can be conditioned on the state of the vertex register but 
do not modify it) and shift (that exchanges the value of the vertex register and 
another register). Nonetheless, we want to restrict the executions of shift to 
(x,y) e E. 

Definition 1. Let 



In the first case we say that Shift succeeds, but in the second case it fails. 

Definition 2 (Problems). Under the restrictions that only static transforma- 
tions and Shift are allowed, consider the following problems: 

— Detect(G): Detect if there is a marked vertex in G; 

— Find(G): Find any marked vertex in G, with the promise that M^fl. 




Detect^ (G) (resp. Detect*-— ^ (G)) will denote the problem Detect(G) with 
the promise that m <G {0, k} (resp. m > k). Similarly, define FlND^(G) (resp. 
FiND ( - fe) (G)) as Find(G) with the promise that m = k (resp. m > k). 

A natural approach to searching on a graph consists in using a random walk. 
Intuitively, a random walk is an alternation of coin flips and shifts. Namely, a coin 
is flipped according to the current state x of the vertex register, its value describes 
a target vertex y, and Shift is a move from x to y. Let p xy be the probability that 
x is shifted to y. Then Shift always succeeds if p xy = whenever (x, y) E. In 
that case, we say that P = (p X y)x.yex is a Markov chain on G. 

We assume from now on that P is an ergodic Markov chain. Therefore P has 
a unique stationary distribution w. We also assume that P is reversible: Tr x p xy = 
T^yPyx, for all x, y G X. To measure the complexity of implementing a random 
walk corresponding to P, we introduce the following black-box operations: 

— Check(M): Check if a given vertex is marked; 

— Setup(P): Draw a sample from the stationary distribution it of P; 

— Update(P): Perform one step of P. 

Each of these black-box operations have the corresponding associated imple- 
mentation cost. We denote by C, S and U the respective complexities of the 
transformations Check(M), Setup(P) and Update(P). 

2.2 Quantum version 

In the quantum case, the problem extends as follows. Let % = C x be a fixed 
Hilbert space with basis (\x)) xe x- Again, a transformation is static if it is con- 
trolled by the vertex register, that is of type J2 xeX \ x )( x \ ® Vx, an d Definition 1 
of Shift is simply extended by linearity. Then the generalization of random 
walks to quantum walks is as follows. 

Definition 3. A quantum walk W on G is a composition of static unitary trans- 
formations and Shift on an invariant subspace of "H^W, the walk space, such 
that Shift always succeeds when W is restricted to its walk space. 

Implicitly we always restrict a quantum walk to its walk space. 

We will only consider quantum walks built from quantum analogues of re- 
versible Markov chains. Thus we extend the operations Check, Setup and Update 
to the quantum setting as follows. Let |0) £ H be a fixed reference state. In the 
following, the first register is the vertex register. 

— Check(M): Map \x}\b) to \x)\b) if x M and \x)\b®l) if x e M, for b = 0,1. 

— Setup(P): Construct the superposition: \ir) = ~^2 xeX ^f^x\ x )'i 

— Update(P): Apply any of V(P), V(P)^ or Shift, where V(P) satisfies 
V(P)\x)\0) = \x)\p x ) = \x) J2 yex y/p^\y), for all x e X. 

Implicitly, we also allow any controlled version of Check(M), Setup(P) and 
Update(P), on which we access via oracle. 

In terms of applications of Shift, Update has complexity 1, and, Setup 
has complexity 0(Sg) (the diameter of G). Nonetheless, in many algorithmic 



applications, the situation is more complex and the number of applications of 
Shift is not the only relevant cost, see for instance [1,2]. 

2.3 Classical hitting time 

From now on we will assume that all the eigenvalues of P are within [0,1]. This 
is without loss of generality by replacing P by (Id + P)/2 if necessary. From the 
ergodicity of P, the eigenvalue 1 has multiplicity 1. The classical hitting time, 
HT(P, M), is defined as the expected number of applications of the Markov 
chain P required to hit a marked vertex when starting from ir. This can be 
used to design a randomized algorithm for Detect and Find based on the 
corresponding random walk. 

Proposition 1. Let k > 1. DETECT^-* 1 ) (G) can be solved with high probability 
and randomized complexity of order 

S + Tx(U + C), where T = max HT(P, M 1 ). 

\M'\=k 

Find(G) can be solved with high probability and expected randomized complexity 
of order 

S + T x (U + C), where T = HT(P, M). 

Let P' be the Markov chain obtained from P by turning all outgoing tran- 
sitions from marked vertices into self-loops. We call P' the absorbing version of 
P. If we arrange the elements of X so that the marked vertices are the last ones, 
matrices P and P' have the following block structure: 

p ._ f Puu Pum\ pi ._ f Puu Pum\ 

\Pmu Pmm J ' \ I J 

where Puu and Pmm are matrices of size (n — m) x (n — m) and m x m, while 
Pum and Pm!7 are matrices of size (n — m) x m and m x (n — m) . 

We first present the matrix characterization of HT(P, M). From the re- 
versibility of P, we get that D(P) — diag( v / 7r)Pdiag( v / 7r)~ 1 is a symmet- 
ric matrix, which is also known as the discriminant of P. The latter was in- 
troduced in [10] for symmetric P, and generalized to reversible P in [12]. 

Set |tt) = ExexV^); Pm - Exgm 71 "* and Pu = T, xe x\M^- The re- 
spective normalized projections of \n) on marked and unmarked vertices are 
\ M ) = J2xeM \[tQvmV) and \U) = J2 xe x\M V^Jpu\ x )- Th en 

BT(P,M)= Pu x (C/|(Id c/ - J D(P) [/[/ )- 1 |;7). 

For simplicity, we will from now on omit the quantity pu in the above defini- 
tion. Indeed, we assume that pu < 1/2, so that the the difference between the 
two expressions is at most a factor of 2. Note that if Pm > 1/2, then there is no 
need for any (classical or quantum) walk to find a marked vertex. 

We now introduce the spectral characterization of HT(P, M). Note that the 
reversibility of P also implies the alternative definition for the discriminant 



(D(P))xy — y/PxyPyx- 



Extending the latter definition of the discriminant to P', we get that D(P') = 
D(P)uu © Idjf. Let \v\), \v 2 ) 7 ■ . ■ , \v n ) be a system of orthonormal eigenvectors 
of D(P') with respective eigenvalues < Ai < A2 < . . . < A n _ TO < A n _ TO+ i = 
. . . = A„ = 1, so that D(P') = J2 k ^k\v k }(v k \- Then one can rewrite HT(P, M) 
as: 

HT(P,M)= J2 TZ^- 

k<n—ni 



2.4 Quantum hitting time 

Quantum walks were successfully used for detecting the presence of marked ver- 
tices quadratically faster than P. Nonetheless, only few is known on the problem 
of finding a marked vertex. Below we illustrate the state of the art. 

Theorem 1 ([10]). Let k > 1. DETECT ( - fc) (G) can be solved with high proba- 
bility and quantum complexity of order 

S + Tx(U + C), where T= max y/HT(P,M'). 

|M'|=fc 

When P is state-transitive and there is a unique marked vertex z (i.e, m = 1), 
HT(P, {z}) is independent of z and one can also find z: 

Theorem 2 ([14,13]). Assume that P is state-transitive. Find (1) (G) can be 
solved with high probability and quantum complexit y of order 
S + T x (U + C), where T = y/WT(P,{z}). 

Using standard techniques, such as in [5], Theorem 2 can be generalized to 
any number of marked vertices, with an extra logarithmic multiplication factor. 
Nonetheless, the complexity does not decrease when the size of M increases, 
contrary to the random walk search algorithm (Proposition 1) and the quantum 
walk detecting algorithm (Theorem 1). 

Corollary 1. Assume that P is state-transitive. Find(G) can be solved with 
high probability and quantum complexity of order 

log(n) x (S + T x (U + C)) , where T = V / HT(P, {zjj, for any z. 

2.5 Szegedy's quantum analogue 

Following the main lines of Szegedy [10], we define a quantum analogue of a 
reversible Markov Chain P. Recall that |0) is an arbitrary reference state in "H, 
and V%P)|a;)|0) = \x)\p x ). Set X = U ® span{|0)} = span{|x)|0) : x e X}, and 
rcf^ = 2J2xex \ x )( x \ ® — Id, the reflection with respect to X. 

Definition 4. The quantum analogue of P is 

W{P) = V{Py ■ Shift • V(P) ■ ref*, 
and its walk space is the smallest invariant subspace that contains X. 



Observe that in [10], the quantum walk is actually defined as (V(P) ■ W(P) ■ 
V(P) 1 ') 2 . Moreover, W(P) requires 3 calls to Update(P), and Shift always suc- 
ceeds when W(P) is restricted to its walk space. 

Szegedy proved the following useful lemma which relates the spectral de- 
composition of W(P) in its walk space to the one of P. Let \vi), \v2), • • • , \v n ) 
be the normalized eigenvectors of D{P) with respective eigenvalues < Ai < 
^2 < • ■ • < A„_i < A„ = 1. From the definition of D(P), observe that |vr) is a 
1-eigenvector of D(P), and therefore we set \v n ) — \n). 

Lemma 1 ([10]L_ Define B k = span{|u fe )|0), W{P)\v k )\0)} , for k ^ n, and 
B n = span{|t!„)|0)}. Then the walk space of W(P) is ($) k B k , where W(P) 
admits the following spectral decomposition: 

- OnB k , k + n: ftf = e ±li ?\ where cos(p k = X k and l ^ +) ^ fc " > = |v fc )|0) 

- OnB n : n n = 1,\&„) = \v n )\0). 

Therefore, we will also call |^„) = |w n )|0) the stationary distribution of W(P). 

3 Finding via Quantum Walk 
3.1 Classical interpolation 

Our starting state is the stationary superposition |i> n )|0) = | vr) 1 0} of W(P). We 
would like to end in its projection onto marked vertices, namely |M)|0), which is 
also the stationary superposition of W(P'). However, in many cases, including 
the 2D-grid, every iteration of W(P') on 1 7r) j 0) may remain far from |M)|0). 

Our approach consists in taking a new random walk, namely an interpolation 
between P and P', This technique is drastically different from the approach 
of [14, 13], and up to our knowledge new. 

Definition 5. The classical interpolation of P and P' is 
P(s) = {l-s)P + sP', 0<s<l. 

This interpolation has some similarities with adiabatic evolutions, where similar 
interpolations are usually defined for Hamiltonians. Here the interpolation is of 
different nature after we take the quantum analogue of P(s). Nonetheless, this 
interpolation still makes sense for adiabatic evolutions, leading to an interesting 
connection between the hitting time and the adiabatic condition [15]. 
Note that P(0) = P, P(l) = P', and P(s) has block structure 

P(s) = ( Puu PuM \ 

[) {(1 - s)P M u (1 - s)P M m + si J " 

Decompose the stationary distribution of P as ir — (iru 7Tm)- Remember that 

Pm — J2xeM n x is the probability to pick a marked vertex from the stationary 

distribution. Then 

*(*) = 1 _ (1 \ s)pu ((1 - s),u n M ) 
is a stationary distribution of P(s), for s e [0, 1]. Moreover one can show that: 
Fact 1 For s £ [0, 1), the Markov chain P(s) is ergodic and reversible. 



For s e [0, 1), let D(s) be the discriminant of P(s). Then D(s) and P(s) are 
similar and therefore have the same eigenvalues. Let \vi(s)), \v2{s)}, ■ . ■ , \v n (s)) 
be a system of orthonormal eigenvectors of D(s) with respective eigenvalues 
< Ai(s) < A 2 (s) < . . . < A„_i(s) < A„(s) = 1: £>(*) - £ fe A fc (s)|u fc (a)X«fc(«)|. 
We also define W(s) = W(P(s)). 

For s = 1, we extend the above for P(l) = f". Recall that |v n (s)) = 
is an eigenvector of D(s) with eigenvalue A„(s) = 1. Observe that \v n (s)) is in 
the two-dimensional subspace spanned by \M) and \U). 

Fact 2 |v n (s)) = cos0(s)|E/) + sin#(s)|M) 7 where 9{s) = arcsin y / ^ s( p 1 " pM) ■ 

Intuitively we want |u n (s)) to have a large overlap on both |{7) and |M), so 
that the algorithm will proceed in two steps: first, map \U) to \v n (s)), using the 
quantum walk (using Update); second, map |t> n (s)) to \M) by projecting onto 
marked vertices (using Check). Therefore, ideally we want s to satisfy sin#(s) = 
cos 9(s) = l/y/2, namely s = s(p M ), where s(p M ) = 1 - jz^j- 



3.2 Quantum circuit for W(s) 

In the following lemma, we assume to know p xx for every x. This is reasonable 
since in practice the probability of self- loops is known. In many cases, it is even 
independent of x. For the rest of the paper, we assume that this is not an obstacle 
(we can assume that one call to Update(P) allows to learn p xx for any x). 

Lemma 2. Assuming that p xx is known for every x, Update(P(s)) can be im- 
plemented with complexity C + U . 

Proof. Since Shift only depends on G, we just need to explain how to implement 
V(P(s)) and its inverse. We simply give an algorithm for V(P(s)). The one for its 
inverse is obtained from the reverse algorithm. From Definition 4, the quantum 
analogue of P(s) is W(s) = F(P(s)) t • Shift • V(P(s) )-rcf x , where V(P(s)) is a 
unitary that maps \x)\0) to \x)\p x (s)) = \x) J2 y ex V 'Pxy{s)\y) ■ We now explain 
how to implement V(P(s j) using one call to V[P) and 2 calls to Check(M). 

Simulation(P, M, s) 

1. Let | £c) 1 0} be the current state 

2. Use a fresh qubit (marking register) in state |0) and call Check(M). 

3. If the marking register is in state |0), call V(P): \x)\p X y}\0} 

4. Otherwise 

(a) Use another fresh qubit in state |0): )a;}|0}|l}|0} 

(b) Apply a rotation of angle arcsin y/s on the fourth register: 
!x)|0)|l)(v / l^|0)+^|l)) 

(c) If the fourth register is |0), apply V(P) on the first two registers, 
Otherwise XOR the first two registers: 

|z> (v 7 ^^} 1 1} |0> + v^l*} 1 1> 1 1}) 

(d) If the second register is \x), apply a rotation of angle 
— arcsin y/s/((l — s)p xx + s) on the fourth register, 

Otherwise do nothing: |a;)|p a! (s)}|l}|0} 

5. Call Check(M) to uncompute the marking register. 



3.3 Hitting time in s 



We define the hitting time in s, which intuitively corresponds to the expected 
time for P(s) to converge to tt(s) from 7r(0). 

Definition 6. HT(,s) = V l( ^ (s)l ^ 2 . 

In particular, note that HT(1) = HT(P, M) is the usual hitting time of P with 
respect to the set of marked vertices M . 

The running time of our quantum search algorithms will depend on HT(s) 
for some particular value s e [0, 1]. To obtain an explicit expression in terms of 
the usual hitting time of P, we show that the maximum of HT(s) is attained at 
s = 1, by providing an explicit expression for HT(s) (proof in Appendix A). 

Theorem 3. HT(s) = sin 4 6(s) • HT(P, M). 



4 Quantum search 

4.1 Algorithm with known parameters 

Theorem 4 (Phase estimation [16,17]). Let W be a unitary operator on 
V. and f£N. There exists a quantum circuit PhaseEstimation(W, t) using 2* 
calls to the controlled operator c—W and 0(i 2 ) additional gates, and acting on 
eigenstates of W as 

2* — 1 
l,m=0 

where e l ' Pk is the eigenvalue of W corresponding to \&k)- 

Observe that a measure of the second register of the output state of 
PhaseEstimation(Wi t)\^k) yields an approximation of the first i-bit of the 
binary decomposition of <^/c/(27r). 

We define our main search algorithm with parameters < p* M < 1 and t £ N. 
Remind that s(p* M ) = 1 — . For suitable parameters, the algorithm outputs 
a marked vertex with high probability, if there is any. 

QuantumWalkSearch(P, M, p*M,t) 

1. Prepare the state | tt) 1 0) 

2. Use a fresh qubit in state |0), apply Check(M) and measure the qubit. 

3. If the outcome is |1), measure the first register (in the vertex basis) and output 
the outcome. 

4. Otherwise, apply PhaseEstimation(H / (s(p^ f )), t) on the registers. 

5. Use a fresh qubit in state |0), apply Check(M) and measure the qubit. 

6. If the outcome is |1), measure the first register (in the vertex basis) and output 
the outcome. 

7. Otherwise, output "No marked vertex" 



Theorem 5. Let p* M ,e\ e [0,1] be such that cos 2 9(s) sin 2 9(s) > e\, where 
s = s(p* M ). Let T > 1 and e 2 € [0, 1] &e such that T > x v /HT(s). ITien, 
QuantumWalkSearch(P, M,}^ , |~logT~|) outputs a marked vertex with prob- 
ability at least e\ — e 2 and complexity of order S + T x (U + C). 
In particular, if \p* M — pm\ < PAf/3 and T > 10-\/HT(P, M), f/ien f/ie success 
probability is at least 1/20. 

Proof. Let i = [log T] . First observe that the complexity analysis is direct since 
QuantumWalkSearch(P, M,p* M ,t) has complexity of order S + 2* x (U + C). 
We now assume that we reach Step 4. Otherwise a marked vertex is already 
found. Then the current state before Step 4 is |t/)|0). Let otk(s) = (U\vk(s)). 
From now on, we omit to write the dependence on s explicitly, when there is no 
ambiguity. 

In Step 4, PhaseEstimation(VF(s), t) is applied on the state 

\U)\0) = a n \v n )\0) + J2 <*k\vk)\0) = a n \* n ) + 4 E ^ (l*fc") + ■ 

k^n k^n 

Theorem 4 shows that PhaseEstimation(T4 A (s), t) maps \H/ n ) to |^„)|0*) and 
maps to |Sf) (£±|0<) + |r,±», where 

;=o m=l J=0 

By definition, (0*1% ) = 0. Then, the probability p to obtain a marked vertex by 
measuring the first register is at least the probability to obtain both a marked 
vertex in the first register and the state |0*) in the last register. Since \Hf n ) = 
\v n )\0), the probability p is therefore lower bounded as 

p > \ an \ 2 \\n M \ Vn )\\ 2 -\Y, Kl 2 (l^ + l 2 + l^l 2 ) = cos 2 sin 2 e £ \a k \Hl 

k^n k^n 

where 77m = SxeM \ x )( x \ * s the projection onto marked vertices and 5 k = \5 k \ = 
\5~k~\- By hypothesis, we already have cos 2 9 sin 2 9 > t\. Therefore, it remains to 
prove that the second term in the RHS is at least — e 2 . 

First, using the definition of 5 k , we get: 5\ = ^ 1^(^/2) - We also 

have by definition of HT(s): 

which together with the above implies that Y^k^n l a *;| 2 ^ 2 — y H ^/^ < e 2- 

We now prove the last part of the theorem. In Appendix B, we prove the 
following fact: 

Fact 3 Let ei < 1/4 be such that 2y/eip M < p M < 2(1 - ^/1~i)pm- Then, 
cos^ 6>(s) sin 2 <9(s) > ei. 

The conditions of Fact 3 are satisfied with e\ = 1/10. Set e 2 = 1/20. Using 
HT(s) < HT(P, M), one can check that the conditions of the theorem are satis- 
fied, and therefore the success probability is at least ei — e 2 = 1/20. □ 



4.2 General case 



At first, assume we have a correct approximation p* M of pm- In that case, even 
if we do not know HT(P, M), we can use the following algorithm, and still find 
a marked vertex with an expected cost 0(y / HT(P, M)). 



QuantumWalkSearch'(P, M,p* M , k) 

1. Let t = 1. 

2. Call k times QuantumWalkSearch(P, M,p* M , t). 

3. If no marked vertex is found, set t 4— t + 1 and go back to step 2. 



Theorem 6. Given p* M such that \p*m — Pm\ < Pm/3, 

QuantumWalkSearch'(P, M,p* M , 28) solves Find(G) wrai/i expected quantum 

complexity of order 

log(P) xS + Px (U + C), w/iere P = ^HT (P,M). 

Proof. Set s = s(p|, ; ). Let t to be the integer such that 7r^/ H ^ s ^ " < 

2 t« < tt^/ 2H ^ (s) , and let P = 2*° = 0( v /HT(s)). By Theorem 5, for 
any t > to, QuantumWalkSearch(P, M,p* M , t) outputs a marked vertex 
with probability at least 1/20. Then, step 3 is reached without finding any 
marked vertex with probability at most p < (1 — 1/20) 28 < 1/4. Moreover, 
QuantumWalkSearch(P, M,p* M ,t) has complexity of order S + 2* x (U + C). 

Let tf be the value oft when QuantumWalkSearch'(P, M,p* M , 28) stops, 
that is, the number of iterations of step 2. Then, the expected complexity of 
QuantumWalkSearch'(P, M,p* M , 28) is of order Nx x S + N 2 x (U + C), where 
iVi is the expectation of tf, and N 2 is the expectation of 2 + 4 + . . . + 2 tf . 

First observe that Ni < t + Efct +iP*~*° = O(t ). For N 2 we get 

to oo CO 

^2 < 2* + ^ ■ 2* = (2 • 2*° - 2) + 2 to J^p* • 2*. 

t=l t=t +l t=l 

Then using the fact that p < 1/4 we finally obtain 

oo 

N 2 < 2 • 2*° + 2*° ^2"* < 3 • 2*°. 

t=i 

This concludes the proof since 2 to = 0(- N /HT(s)) and HT(s) < HT(P, M). □ 

For the general case we get two possible situations, depending whether a 
lower bound p m - m on pu and/or an upper bound HT max on HT(P, M) is given. In 
particular, for FiND(G)'- fc \ we can set p min = min M /.| M /| =fe p M i and HT max = 
max M , : | M ,| =fe HT(P,M'). 

Theorem 7. Given p m [ n < pu, Find(G) can be solved with expected quantum 

comp lexity of ord er 

v/log(l/p min ) x [log(T) x S + T x (U + C)] , where T=y/WT(P,M). 
Moreover, i/HT max > HT(P, M) is also given, then Find(G) can be solved with 

expected qu antum comp lexity of order 

v/log(l/p min ) x [S + P x (U + C)] , where T = yJWT mm . 



Proof. We simply prove the first part of the theorem. The second 
one is similar using QuantumWalkSearch(P, M,p* M , T) instead of 
QuantumWalkSearch'(P, M,p* M , 28). 

From Theorem 6, it is enough to have a good approximation p* M oipM, such 
that ip* M /A < pm < 3p* M /2. Moreover, since p m ; n < pu < 1/2, this condition 
will be satisfied for some p* M e {(2/3) x 2~ l : I = 1, . . . , [log(l/p m i n )J }• 

Let us incorporate step 2 of QuantumWalkSearch'(P, M,p* M , 28) into a 
loop on the [log(l/p m ; n )J possible values of p* M . Then the analysis is basically 
the same, except that now the complexity of step 2 is multiplied by a factor of 
order log(l/p m ; n ). This loop is then implemented a la Grover using the approach 
of [18], reducing the multiplication factor to ^/log(l/p min ). □ 

Theorem 8. Given HT max > HT(P, M), Find(G) can be solved with expected 

quantum complexity of order 

log(l/p M ) x [S + T x (U + Q] , where T = v / HT max . 

Proof. We now use QuantumWalkSearch(P, M,p* M , t) with t = 
[log VHT max ] , and perform a dichotomic search for an appropriate value 
of p* M . This dichotomic search uses backtracking since the branching in the 
dichotomy is with bounded error. This is described in [19]. 

Initially we set a — and b = 1. Then for testing the cur- 
rent value of p* M = (a + b)/2, we run a constant number of times 
QuantumWalkSearch(P, M 7 p* M , t). If a marked vertex is found we stop. Oth- 
erwise, if PhaseEstimation(VF(s(p^ // )), t) outputs a minority of Os, we set 
a = p* M , otherwise we set b — p* M . □ 

4.3 Application to the 2D-grid 

We consider in this section only the complexity in terms of the number of uses 
of Check and Shift. Consider a random walk on the 2D-grid of size y/n x t/n, 
with self-loops. The previous best known quantum complexity of FiND(G)( fc ) 
and FiND(G)(- fe ) was 0( v / ^(logn) 3 / 2 ), from Corollary 1. Since the grid is a 5- 
rcgular graphs (4 directions and 1 self- loop), P is symmetric, and therefore the 
stationary distribution of P is uniform, and we simply have pu = m/n. Then 
Setup is realized with y/n uses of Shift, and HT(P, {z}) — 8(nlogn), for any 
z. Therefore we get the following corollary of Theorem 5 and Theorem 7, by 
upper bounding HT(P, M) — O(nlogn). 

Corollary 2. Let G be the 2D-grid of size y/n x y/n, and let k > 1. Then 
FiND(G)( fe ) can be solved with expected quantum complexity 0(\/n logn), and 
FlND(G)(- fe ) with expected quantum complexity 0{\Jn x logn x log(n/fc)). 
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A Proof of Theorem 3 



Let us first consider the derivative of HT(s). 

Lemma 3. The derivative ofWI(s) is related to HT(s) as: 

ds l-s(l-p M ) 

Note that Theorem 3 follows directly from this lemma, as we may obtain an 
explicit expression for HT(s) by solving the differential equation for HT(s) and 
using HT(1) = HT(P, M) as boundary condition, which yields 

HT(s) = ?M j HT ( P > M ) = sin4 6 ( s ) ■ HT ( P ' M )- 

(1 - s(l -p M )) 

Therefore, it remains to prove Lemma 3. 

In this section, we will often use f(s) as a shorthand form of j^f(s)- Let 
II m = J2xgm \ x )( x \ De the projector onto the m-dimensional subspace spanned 
by marked vertices, and let {X, Y} := XY + YX be the anticommutator of X 
and Y. 

Before proving Lemma 3, let us consider the following claim. 
Proposition 2. D(s) = —L^{n M ,I-D(s)}. 



Proof. Using the definition D(s) = y/P{s) o P T (s) and the block structure of 
P(s), we have: 



D(s) 



^PUUOP^U ^(l-s){P UM oPl 
J (I - s)(P MU o PJ m) (1 _ s )^p MM o P T M + si , 



and hence 



D(s) = 



PUM O P M U 



2^n^s 

J^= s \/PmU°Pum 1 ~ V ' P MM O P MM 

Observe that D(s) + 2 (i_ s ) {^m, -P( s )} = jhr^ Hm, which implies the claim. □ 
Define 

Then from Fact 2 and the observation that A(s)\v n (s)) — 0, we get 
Fact 4 sm9(s)A(s)\M) = - cos 0(s)A(s)\U). 
We are now ready to prove Lemma 3. 



Proof (Proof of Lemma 3). In this proof we will often omit to write the depen- 
dence on s explicitly. From definitions of HT(s) and A, we have 

HT = (U\A\U). 

Note that A = B^ 1 - LJ n , where 

B = ld-D + n n , n n = \v n )(v n \. 

For any invertible matrix M we have Jj(M _1 ) = — M^MM^ 1 . Therefore, 

A = -B- 1 (-b + ti n )B- 1 -tl n . 

Hence, we have ^ HT = hi + h 2 + /13, where 

hi = (UIB^DB^IU), 

h 2 = -(u\B- 1 n n B- 1 \u), 
h 3 = ~{u\n n \u). 

Let us evaluate each of these terms separately. We can use Proposition 2 and 
the definition of B to express the first term hi as follows: 

2(l-s)hi = {U\B- 1 {n M ,I-D}B- 1 \U) 

= {u\B- 1 ({n M ,B}-{n M ,n n })B- 1 \u) 

= (u\{B-\n M }\u) - (uib-^ltm^^b-^u). 

Note that LTm\U) = 0, thus the first term vanishes. Also note that B _1 |w„) = 
\v n ) and B~ x is Hermitian, thus 

2(1 - s )hi = -2{u\B- x n M n n \u). 

Note from Fact 2 that LTM\v n ) = sin#|M). Moreover, using Fact 4 we can 
simplify equation (A) even further: 

2(1 - s)hi = -2cos9sm9(U\B- 1 \M) 

= -2cos0sin0(t/|(,4 + i7„)|M) 

= -2cos<9sin<9(J7|,4|M) - 2 cos 2 9 sin 2 9 

= 2 cos 2 9((U\A\U) - sin 2 9). 

Let us now consider the second term h 2 : 

h 2 = -(U\B- 1 il n B- 1 \U) 

= -cos9((v n \B- 1 \U) + (UIB- 1 ^)) 

= -2cos9(U\(A + n n )\v n ) 

= -29 cos 9 (-sin 9 (U\A\U) + cos 9 (U\A\M)) 

/ cos 2 9 

= -29cos9 l-smO(U\A\U) - —^(U\A\U) 



where we have used the fact that 

\v n ) = 0(-sin0|Z7) +cos0|M)), 
and II n \v n ) = 0. Similarly, for the last term I13 we have 
h 3 = (U\fl n \U) = 26 cos 6 sin 0. 
Putting all the terms back together, we have 

d = ^l,, u , A , ir ._^ n 2 ^ +2 e^l(u\A\U) + 26cos6sm8 
ds 1 — s sin 6* 

cos6> / • cos#sin0\ , . A . T , „ . . /„; cos0sin#\ 

= -r-^ 2(9 + — f/L4f/ +cos6»sin6» 2(9 

sin \ 1 - s J \ 1 - s J 

From the definition of 8, we directly get 



2 ^ _ cos 9 sin fl _ \/pm(1 ~Pm) cosfl _ / (1 - s)(l - Pm) 

l-s y/1 - s(l - s(l-p M )) ' si 116 * V Pm 

which together with the above implies the lemma. □ 
B Proof of Fact 3 

Since s — 1 — , by definition of 6(s) we have 

sing( S )= ( 1 ;Pm)pm and cos0(s)= Pm(I-Pm) 



p M + (1 - 2p M )p* M ' Pm + (1 - 2p M )p* M ' 

From the above expression, we have sin6>(s) > y 7 ^! if and only if 

" ViT+(l-2V?T)pM- 
Since ei < 1/4 and pu < 1/2, the denominator is upper bounded as 

Therefore, p* M < 2(1 — y/ei)pM implies sin f9(s) > -yje[. 
Furthermore, we also have cos#(s) > ^/ej if and only if 



PM > ^ 



(i-^r 1 )-a-2^r 1 )p M ' 

where the denominator is lower bounded as 

1 - 2^1 _ 1 



(1 - v^T) - (i - 2V^)pm > (1 - v^T) - 

Therefore, p* M > 2 s fe{p M implies cos#(s) > y/ei. 
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